Apparent hysteresis in a driven system with self-organized drag 
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Interaction between extended defects and impurities lies at the heart of many 



physical phenomena in materials science. Here we revisit the ubiquitous prob- 
lem of the driven motion of an extended defect in a field of mobile impurities, 
which self-organize to cause drag on the defect. Under a wide range of exter- 
nal conditions (e.g. drive), the defect undergoes a transition from slow to fast 
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O ■ motion. This transition is commonly hysteretic: the defect either moves slow 



or fast, depending on the initial condition. We explore such hysteresis via 



a kinetic Monte Carlo spin simulation combined with computational coarse- 
graining. Obtaining bifurcation diagrams (stable and unstable branches), 
we map behavior regimes in parameter space. Estimating fast-slow switching 
times, we determine whether a simulation or experiment will exhibit hysteresis 
depending on observation conditions. We believe our approach is applicable 
to quantifying hysteresis in a wide range of physical contexts. 
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Driving an extended defect or domain wall in a system that contains stationary or mobile 
impurities with which the domain wall interacts lies at the heart of many physical phenom- 
ena, including the motion of grain boundaries in polycrystalline materials [1], the motion 
of vortex lines in dirty superconductors [2,3], driven charge-density waves [4], the motion 
of dislocations in impure metals undergoing plastic deformation (i.e., Portevin-LeChatelier 
effect [5]), ferroelectric domain wall dynamics [6], and stick-slip phenomena in tribology [7]. 
The ability to quantitatively model such systems holds the key to designing a wide range 
of devices and optimizing their operating conditions. While understanding the transition 
between the slow and fast kinetic regimes is vital, the dynamics become particularly difficult 
to characterize in the neighborhood of such a transition (neutral stability). 

When equations that describe the macroscopic behavior of a system are known analyti- 
cally, an established set of mathematical/computational tools is available for analyzing the 
dynamics (e.g., locating and characterizing transitions). Many "real world" applications 
are too complex to characterize in terms of simple equations. Microscopic simulations are 
becoming an increasingly common approach to analyze situations where such macroscopic 
models are not available analytically. This is the case for driving a domain wall through a 
field of mobile impurities: although analytical models exist, they do not accurately and/or 
fully predict the range of observed behavior. In such systems, there is a "specter" of hys- 
teresis in the transition from slow to fast motion - sometimes it is seen [8] and sometimes 
not [9] . Classical theories of such phenomena [1] are usually one-dimensional and show that 
the domain wall velocity vs. driving force curve can be multi- valued, i.e., suggest the exis- 
tence of hysteresis. In this brief report, we examine this phenomenology - try to understand 
its roots, and argue that it is both the physics and the observer that determines whether 
hysteresis will appear in experiments or simulations. 

We illustrate the issues arising in considering hysteresis within the context of the classical 
solute drag effect, in which diffusing impurities interact with a migrating domain wall and 
impede its motion. In particular, we employ a variant of the Ising model in which the 
domain wall moves under the influence of an external field, and interstitial impurities diffuse 
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and are attracted to the domain wall [9]. The energy of the model system is 7i/k B T = 
—J/2 SiSj — hJ^Si — E /4:J2 e a\J2 s j\, where Sj = 1 in one domain and —1 in the other, e a 
is zero or one depending on whether an impurity occupies interstitial site a or not, E > 
implies that the impurities are attracted to the interface, and J > and h scale the domain 
wall energy and the driving force, respectively. The first term represents the total energy 
of domain walls in the pure system, the second term accounts for the external drive, and 
the last term describes the interaction between the impurities and the domain wall. The 
parameters chosen for the present study were J = 2.0 and E = [4.0, 7.0] and the bulk 
impurity concentration was set to 1% (cj mp = 0.01). 

Initially, the 2-d computational domain contains a single, straight domain wall and ran- 
domly placed interstitial impurities. The system evolves by randomly choosing spins to 
flip or impurities to exchange with neighboring empty sites, with probability proportional 
to their relative density and the impurity diffusivity m, as described by Mendelev, et al. 
[9]. An attempted change is accepted if it lowers the energy of the system (AH < 0) or, if 
AH > 0, it is accepted provided that e _AH//fcsT < R, where R is a random number uniformly 
distributed between and 1. The simulation clock advances by one unit when all of the 
spins have been sampled once, on average. When the domain wall gets close to one end of 
the simulation cell, more sites are added to that end and removed from the opposite end, 
thereby making the length of the system L effectively infinite [9]. The results presented 
below are for a relatively small simulation cell width (W = 32). This was done to prevent 
excessive domain wall roughening, which tends to smear the slow-fast domain wall transi- 
tion. The effect of system size is very important for the observation of hysteresis, but will 
be the subject of a future report. 

The domain wall velocity results from a balance between the external field, which controls 
the propagation velocity in the pure system, and the diffusing impurities, which locally 
impede domain wall migration. In the limit of large drive and weak domain wall-impurity 
interactions, the impurities have little influence on the domain wall motion. Conversely, 
decreasing the drive and increasing the attraction of the impurities to the domain wall E 
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leads to very slow domain wall dynamics, due to the strong segregation of the impurities 
onto the domain wall. However, as will be demonstrated below, the drive, impurity-domain 
wall interaction, impurity mobility and the average impurity concentration may conspire 
to give rise to an intermediate regime: a regime where the domain wall is captured by the 
impurities, then escapes and travels nearly unimpeded until it is recaptured. This is the 
"jerky" (stick-slip-like) domain wall motion regime, where hysteresis may be observed. 

Figure 1 shows the mean domain wall position R vs. time £ for h — 0.12 and several values 
of the domain wall/impurity interaction strength E . These curves become qualitatively 
different as E is changed. In particular, for E = 4.0, the domain wall propagates smoothly, 
while increasing E from E = 5.5 to E = 6.5 produces "jerky" motion, in which the domain 
wall slows dramatically for long time intervals before breaking free again. Examination of 
Fig. 1 also shows that intervals of slow motion are strongly correlated with the presence of 
a large number of impurities on the domain wall, N, and vice versa. Transitions in domain 
wall speed are strongly correlated with a jump in the number of impurities on the domain 
wall. 

It is instructive to consider how the changes in R(t) are reflected in the distribution 
of the number of impurities on the domain wall P(N), shown in Fig. 2, obtained from 
long-time simulations at several Eq. For relatively weak domain wall/impurity interactions 
(E = 5.25), P(N) is unimodal - the system is (almost always) in the fast regime. As the 
interaction strength is increased, P(N) becomes bimodal as the system transits between the 
slow and fast regimes (i.e., jerky motion). Finally, for sufficiently strong interactions, P(N) 
is again unimodal, with the system effectively trapped in the slow state. 

In the "jerky" regime, observations of the system can be characterized by two distinct 
time scales, namely ts, the characteristic time that the domain wall motion remains slow 
(effectively trapped, many impurities on the domain wall) before it breaks free; and tp, the 
characteristic time that the domain wall motion is fast (few impurities on the domain wall) 
before again being captured by the impurities. If our observation time to is smaller than 
both these characteristic times, the system will appear to us in a typical simulation to be 
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either only fast or only slow, depending on the initial condition (suggesting hysteresis). On 
the other hand, if t Q is much greater than these times, the domain wall appears to move at 
a constant speed that is an appropriately weighted average of the fast and slow speeds. The 
motion appears to be "jerky" only if to is in the range of the two characteristic times. 

In what follows, we focus on the variable N(t) and demonstrate that a quantitative 
description of the system can be developed based upon its dynamics. Assuming that N(t) is 
an effective "slow variable" for the system, standard arguments based on the Mori-Zwanzig 
projection operator formalism (see, e.g., Ref. [10]) lead to a stochastic description of the 
system in terms of a Langevin equation for N(t). Alternatively, we may consider a one- 
dimensional Fokker-Planck (F-P) equation, describing the statistics of the stochastic motion 
of a "particle" in a two-well potential; here the "particle" corresponds to the instantaneous 
number of impurities on the domain wall N(t), and the effects of all other degrees of freedom 
are incorporated into the effective potential. The effective F-P equation for the corresponding 
probability distribution P(N, t) is [11] 

dP(N,t) 



dt 



P(N,t), (1) 



dN y ' dN 2 

where V(N) = \im At ^ (AN) /At, 2D(N) = lim At ^ ([AA^] 2 )/At, and AN = N(t + At) - 
N(t). At stationarity, P(N) oc exp(— $(A^)//csT), where $(iV) denotes the effective poten- 
tial and is determined from V and D as 

$(N) f N ,„,V(N' 

' = const. 1 ' v ' 



dN 



k B T Jo D(N' 

Note that, in the limit D = const., Eq. (1) becomes equivalent to the Langevin equation 
:>t = ~ fc^fzf 7- '"^' wnere V 1S a random, stochastic variable that is Gaussian distributed with 
mean (n) = and variance (n(t)n(s)) = 2D5(t — s). Most importantly, our MC simulations 
can be employed to extract both V(N) and D(N) "on demand" as follows: choose a value 
of No, locate several instances when it appears in a long time simulation, tabulate the 
subsequent values of N within a fixed time interval At (here At = 600), and then average 
over these segments and estimate the rate of change in the mean (V) and the variance (2D) 



for the number of impurities on the domain wall. This is repeated for a grid of N values 
sufficient to numerically evaluate the integral in Eq. (2) to the desired accuracy. Similar 
methods for estimating V and 2D from a stochastic time series has been presented in Refs. 
[12,13]. 

The result of this analysis is shown in Fig. 3, where we plot &(N)/k B T for E = 6.0 along 
with the potential obtained by directly constructing the probability distribution P(N) from 
the time series and employing P(N) ~ exp(&/k B T). Indeed, §(N)/k B T has a double-well 
structure, with minima at N B ~ 4 and Ns ~ 24, and an (unstable) saddle transition point at 
Nu = 11. Furthermore, the effective potential extracted from the F-P analysis shows good 
agreement with that obtained directly from the time series for N > 6. The discrepancy 
observed for small N can be rationalized by arguing that $(JV) is no longer effectively one- 
dimensional. In this regime, additional degrees of freedom, such as the domain wall shape, 
are no longer slaved to the single variable N(t). 

The form of the effective potential $(iV) suggests that there is a single effective barrier 
or saddle that must be overcome for the system to move from the slow (pinned) state to 
the fast one (or vice versa). The waiting time between transitions between states depends 
on the shape of $(JV) and the kinetic coefficient D(N). Following Kramers [14], we can 
estimate this time r, as 

T ~ 2ixk B T e A<f/fc B T / 3 n 

~ D^"(N min )\<Z>"(N saddle )\ 

where 2D = D(N min ) + D(N sa ddie)- For the data represented in Fig. 3 corresponding to 
E = 6.0, the average time needed for a transition from fast to slow and vice versa is 
t m 4 x 10 4 and r ps 4 x 10 5 , respectively. These are each within a factor of 3 of the waiting 
times estimated from a long MC simulation. 

Armed with $(iV) for several values of E Q , we now construct the domain wall speed 
vs. impurity interaction strength, E , bifurcation diagram, as follows. We determine the 
average domain wall speed as a function of N, V = V(N), for each E from our numerical 
simulations. Then, we determine the velocity on the slow branch as Vs = V(N = Ns), on 
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the fast branch as Vf — V(N — N F ), and on the unstable branch as Vu — V(N — Nu), 
where N s , N F , and Ny) from plots such as those in Figs. 2 and 3 for each value of E . The 
results are shown in Fig. 4, where we plot the stationary domain wall velocity vs. interaction 
strength bifurcation diagram, including multiple stable branches and the unstable branch. 
Alternatively, we can search for the stationary values of A" for which dN/dt = using 
standard root-finding methods. dN/dt can be estimated as N ( t -°)- N ( t - t ) by performing a 
series of short simulations initialized with N(t = 0) impurities on the boundary for a time 
? [13]. 

To summarize, we obtain both the stable and unstable branches of the bifurcation di- 
agram, as well as an estimate of the characteristic time r over which the domain wall will 
switch from fast to slow or vice versa. Therefore, hysteresis can be expected for obser- 
vation times t <C r, and the bifurcation diagram becomes single-valued everywhere for 
observation times to 3> t. In the latter limit, the observed domain wall velocity will be 
V(E ) = Jo°° dNP(N;E )V(N;E ), as shown in Fig. 4. Increasing either the impurity 
difiusivity or the heat of segregation (or both) leads to more effective pinning, and thus 
increasingly jerky domain wall dynamics. In Fig. 4 we show a morphological diagram for 
the propagation mode of the domain wall as a function of m and E for t Q = 10 4 . 

Combining a coarse-grained computational approach with an effective Langevin/Fokker- 
Planck description for the stochastic dynamics of the impurity density along the domain wall 
provides a useful and computationally efficient method for determining the apparent station- 
ary states (both stable and unstable) of the system, as well as for computing the transition 
rates between the stable states. The transition rates, together with the apparent stationary 
states, provide a clear picture of the dynamics of the coupled domain wall-impurity system 
and conditions for determining when hysteretic behavior will be observed. We believe that 
this approach can be used to quantifying the emergence of hysteresis in a wide variety of 
physical systems. Whether hysteresis will be observed in a particular experiment or simu- 
lation depends, of course, on the physics. However, it does not only depend on the physics; 
it depends on the time scale of the observation and on the (length, ensemble) size of the 
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system. In this sense, hysteresis is in the eye of the beholder. 
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FIGURES 

FIG. 1. (a) Position of the domain wall as a function of time for several values of the heat 
of segregation, Eq. (b) The number of impurities N(t) on the domain wall, corresponding to the 
Eq = 6.0 simulation in (a). (c,d) Configurations of the domain wall and impurity positions at the 
times marked "X" and "O" in (a), respectively. 

FIG. 2. Steady-state probability distributions P(N) as a function of Eq, determined from long 
MC runs. 

FIG. 3. The effective potential <3>(iV) from Eq. (2) and directly from the long-time simulation 
data of Fig. 2 with E = 6.0. 

FIG. 4. (a) The effective domain wall velocity versus Eq bifurcation diagram at m = 0.25. The 
solid and dashed lines are guides to the eye. (b) A diagram of the different regimes of behavior 
as a function of impurity mobility m and the impurity-domain wall interaction strength Eq; blue 
denotes smooth, gray jerky, and red effectively pinned propagation. Symbols are from measured 
data with "F" corresponding to fast, "J" to jerky, and "S" to slow (pinned) behavior, respectively. 
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